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The Green's function formalism for neutrino mixing is presented and the exact 
oscillation formula is obtained. The usual Pontecorvo formula is recovered in the 
relativistic limit. 



1 Introduction 

We report on,the Greenis-function formalism for neutrino field mixing recently 
presented in 111 (see alsoS'0). The result is an oscillation formula which differs 
from the usual oneta in the non-relativistic region. We get, together with the 
"squeezing" factor of the amplitude found in ref. 13, also an additional term with 
a different oscillatory frequency. This last feature is particularly important 
since it shows that resonance is possible also in vacuum for particular values 
of the masses or of the momentum, thus leading to a suppression or to an 
enhancement of the conversion probability. 

We consider two Dirac neutrino fields v e and (space-time dependence 
suppressed). The "flavor mixing" transformations are 

v e (x) = cos9 + vi(x) sinf? 

— —vi(x) sva9 + V2{x) cosf? , (1) 

where 9 is the mixing angle. v\ and 1/2 are explicitly given by 

Vi(x) = ±= [^e-^'Xi + V-Ki^'^i] e lk ' x , i = 1, 2 . (2) 

We use t = Xo, when no misunderstanding arises. The vacuum for the ai and 
Pi operators is denoted by |0)i,2: JO) 12 = /3£ J0)i2 = 0. The anticommu- 
tation relations, the completeness and orthonormality relations are the usual 



i 



ones. In order to cimiimvent the difficulty of the construction of a Fock space 
for the mixed fields 130, it is useful to expand the flavor fields v e and in the 
same basis as v\ and 1^2; e.g. v e is given by 

u^{x) = G~ l {e,t) v^{x)G{9 1 t) 

^ k,r 



where G(9,t) = exp ^9 Jd 3 x (jA{x)v2{x) — v \(x)vi(x)^j 
the mixing transformations (111)13. The flavor annihilation 



is the generator of 

the mixing transformations The flavor annihilation operators can be now 

explicitly given. For example, the electron neutrino annihilator is 

o£ e (i)=cos# o£ ]+ 
sinfl £ (ulW^e-^-^ay + u^vt^e^ + ^(3%) . (4) 

Notice that_it has contributions from ct\, cti but also from the anti-particle 
operator @\ □ since spinor wave functions for different masses are not orthogo- 
nal. In the more traditional treatment of mixing, the @\ contribution is missed 
since the non-orthogonality of the spinor wave functions is not considered. 

We can show that when the two point Green's function for the mixed 
fields v e , are constructed by using the vacuum |0 >i.2, then the "survival" 
probability amplitude, say, of an electronic neutrino state in the limit t — > + 
is computed to be V e e(k,0 + ) — cos 2 9 + sin 2 9 |t/k| 2 < 1, which is clearly not 
acceptable since, of course, it should be lim t _, + Veeif) = 1- 

Here \U^\ 2 is calculated from the spinor basis and ita explicit form is given 
inEl. For different masses and k ^ 0, |[/k| is always < 10, and we will also use 
|Vk| = v/l — |?7k| 2 . |C^k| 2 — * 1 in the relativistic limit k 3> ^/mlmj. 

The above contradiction shows that the choice of the state \Q)i,2 in the 
computation of the Wightman function is not the correct one. The problem 
is in the fact that the transformation (|]) does not leave invariant the vacuum 
1 0) 1.2 • The mixing generator induces on it a SU(2) coherent state structure, re- 
sulting in a new state, \0(9, t)) e _^ = G^ 1 (9, 4)10)^2 , which is the flavor vacuum 
for the flavor operators a e/fl , f3 e/ jB: a 1 k e/fi (t)\0(t)) etf _ l = P^ e/fJi (t)\0(t)) e ^ = 0. 
An important feature of the flavor vacuum |0) eiM (and of the relative Fock 
space) is its non-perturbative nature, resulting in the unitary inequivalence 
with the "perturbative" vacuum 1 0) 1^2 , in the infinite volume limit □. Notice 
that the squared modulus of the survival probability amplitude reproduces the 
Pontecorvo oscillation formula in the relativistic limit. 

We show below that the correct definition of the Green's functions is the 
one which involves the non-perturbative vacuum \0) e ,fj,- 
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2 Green's functions for flavor neutrinos 



In the case of v e — > f e propagation, the relevant Wightman function is (we 
use xo =t,y = 0) iG>^(t,x;0,y) = ei/1 (0|i/°(i,x) ff (0, y)|0) e , M . It can be 
conveniently expressed in terms of anticommutators at different times as 

iG>^(k,t) = 

r 

Here a^ e stands for a k * e (0) . The corresponding transition amplitude is 

P e r e (k,i) = cos 2 6> + sin 2 6> r|£/ k | 2 e -*(«*.!»-««..i)* + |v k |2 e <(««..a+"*,i)tl _ ( 6 ) 

We thus find that the probability amplitude is now correctly normalized: 
lim t ^ +Vee(k,t) = 1, and one can show that V ee , V^, "P^ e go to zero in the 
same limit t — > + . Moreover, 

|P e r e (k,t)| 2 + |7* (k,t)| 2 + |7^(k,t)| 2 + (k,t)| 2 = 1 , (7) 

as the conservation of the total probability requires. Notice that in the per- 
turbative case, there were only two non-zero amplitudes, i.e. V ee and V^. 

At time t = the one electronic neutrino state is (momentum and spin 
indices dropped) \v e ) = a£|0) e>M . In this state a multiparticle component is 
present, disappearing in the relativistic limit k 3> ^Jm\ra2 , where the Pon- 
tecorvo state is recovered. Its time evolution is given by \v e (t)) = e~ tHt \v e ) 
and in the flavor basis this state is found to be 

\v e {t)) = h(t) 4 + m (t) a t + m (t) aWA + m{t) ajaj./^] |0) e , M • (8) 

Here the T)(t) are coefficient satisfying the normalization condition |??i(i)| 2 + 

\m{t)\ 2 + \m{t)\ 2 + \m{t)\ 2 = 1. 

Notice that |0) e)A , is not eigenstate of the free Hamiltonian H; it "rotates" 
under the action of the time evolution generator: \0(t)) e ^ = e~ lHl - 2t |0) e/i . In 
fact one finds limy^oo e ,/x(0 | 0(i)) ej([i = 0. Thus at different times we have 
unitarily inequivalent flavor vacua (in the limit V — > oo): this is not surprising 
since it is direct consequence of the fact that the flavor states are not mass 
eigenstates and therefore the Poincare structure of the flavor vacuum is lacking. 

Finally, the charge operators are Q e / M = a \/n 0le /i JI - ~ Pe/pPe/v We have 
e,n(0{t)\Q e / M|0(*))e l( u = and charge conservation is ensured at any time: 
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(v e (t)\{Q e + Q^) \v e (t)) = 1. The oscillation formula for the flavor charges 
then readily follows 

P„.^„.(k,i) = 



sin 2 (20) 



|f/k|2 s . n2 / A + |14|2 sin2 ( f 



sin 2 (20) 



2 '^).|^ 8in 2^fe,2+^,l 



It/kl 2 sin 2 



(9) 



(10) 



This result is exact and includes the previous result of momentum dependent 
oscillation amplitude of refsH. Notice that the additional contribution to the 
usual oscillation formula, does oscillate with a frequency which is the sum of 
the frequencies of the mass components. In the relativistic limit k 3> yjra\rri2 
the traditional oscillation formula is recovered. 
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